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Abstract. Since its launch INTEGRAL has devoted a significant fraction of its observing time to the Galactic 
plane region. We present the results of INTEGRAL observations of the inner spiral arms of the Galaxy (Norma, 
Scutum and Sagittarius) with attention to high mass X-ray binaries. The increased sensitivity of the survey and 
possibility to detect absorbed sources significantly enlarged our sample of HMXBs in comparison with previous 
studies. For some absorbed sources we present for the first time high quality broadband (1-100 keV) energy spectra. 
We confirm pulsations of the X-ray flux of IGR J16358-4726, discovered by CHANDRA, and report discovery of 
pulsations with a period of 228 ± 6 sec from IGR J 16465-4507. We confirm that the Galactic high mass X-ray 
binary distribution is significantly different from the distribution of low mass X-ray binaries concentrated in the 
Galactic Center. A large fraction of detected high mass X-ray binaries are absorbed sources. 
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1. Introduction 

One of the most interesting results obtained by the 
INTEGRAL observatory is the discovery of a large 
population of highly photoabsorbed Galactic X-ray 
sources. The first source discovered with INTEGRAL 



(IGR J16318-4848, lUourvoisier et al. 2003|> was shown 
to have a large intrinsic absorption in the X- 



ray spectrum (e.g. 


Murakami, Dotani, & Wijnands 2003 


Ide Plaa et al. 20031 


Matt L Guainazzi, 20031 Revnivtsev 



et al. 2003a). After this source several more were dis- 
covered (e.g. IGR J16320-4751=AX J1631. 9-4752, IGR 
J16358-4726, etc.) which also have high photoabsorption, 
strongly exceeding that of the interstellar medium. It was 
proposed that these sources can be members of a previ- 
ously poorly known population of Galactic X-ray sources 
with high intrinsic absorption IjRevnivtsev et al. 2003a|l . 
Based on the photometric observations of the source IGR 
J16318-4848 these authors also suggested that the source 
could be a binary with a giant or supergiant companion 
with a strong stellar wind, that was later confirmed by 
precise spectroscopic observations of the optical counter- 
part HFilliatre fc (Jhaty 2004| ). 

The analysis of properties of newly discovered ab- 
sorbed sources showed that they are likely high mass 
X-ray binaries (HMXBs) with early type companion 
stars, and some of them contain accreting X-ray pulsars 
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ilPatel et al. 20041 IWalter 20041 ILutovinov et al. 2005afl . 
Remarkably these sources are confined in a relatively small 
region of the sky, close to the tangent to the Norma 
Galactic spiral arm. Based on a previous study of the pop- 
ulation of high mass X-ray binaries in our Galaxy (Grimm 
et al. 2002) it was proposed (Revnivtsev 2003) that such a 
concentration of absorbed HMXBs is the result of an en- 
hanced star formation rate in this region and is an obser- 
vational appearance of the Galactic spiral arm structure. 

A similar enhanced population of HMXBs (more 
specifically, accreting X-ray pulsars) in the regions of spi- 
ral arms was emphasized by Koyama et al. (1990). Basing 
on the GINGA observations it was shown that in the re- 
gion of the Scutum spiral arm tangent there exists a large 
number of accreting X-ray pulsars, likely variable, that 
was considered as an indication that they are in binary 
systems with B[e] companion stars. A subsequent study 
of these systems with the ASCA observatory showed that 
many of them are absorbed ( |Sugizaki et al. 2001) ) , and the 
comparison of the measured value Nh with the interstel- 
lar Galactic absorption indicates that the observed X-ray 
photoabsorption is likely intrinsic. 

Previous studies of Galactic X-ray sources (e.g. 
Koyama et al. 1990, Negueruela 2004) suggest that sources 
with significant intrinsic photoabsorption in their spectra 
can be unambiguously classified as HMXBs, while the ab- 
sence of large photoabsorption does not allow a reliable 
classification. 
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Such a difference in observational appearances of high 
mass and low mass X-ray binaries can be understood if we 
remember that in the majority of high mass X-ray binaries 
the compact objects accrete matter from the stellar wind 
of a young companion star, and therefore the dense stellar 
wind can create an envelope which will result in signifi- 
cant photoabsorption within the binary system. The best 
known examples of high mass X-ray binaries in our Galaxy 
are Cyg X-3, Vela X-l, GX301-2 and 4U1700-37. They all 
show the strong photoabsorption in their X-ray spectra. 

In the case of low mass X-ray binaries (LMXBs) the 
main accretion goes through the inner Lagrangian point 
and therefore practically no obscuring matter is available 
in such systems (except for the rare class of dippers or ac- 
cretion disk corona sources, see e.g. Kahn 1982). From this 
reasoning we can expect that any sample of HMXBs se- 
lected from soft X-ray data might lack systems with high 
intrinsic photoabsorption and therefore will lack a con- 
siderable number of HMXBs. In order to construct more 
uniform sample of HMXBs a hard X-ray (> 10 - 20 keV) 
energy band survey is important. 

The international gamma ray observatory INTEGRAL 
( Wink ler et al. 2003"jl is useful for the search of absorbed 
sources. The telesc ope IBIS/ISGRI l|Ubertini et al. 20031 
ILebrun et al. 2003| has a large effective area (~1000 
sq.cm.) and works in the hard X-ray energy range (> 20 
keV), that largely excludes any influence of the photoab- 
sorption of X-ray sources on the detection efficiency. It has 
a large field of view (~ 29° x 29°) and good localization 
accuracy (~ 2 — 3' ). 

The first two years of INTEGRAL observations 
provided more than 20 million seconds of data. The 
deepest hard X-ray (~20-60 keV) surveys to date of 
the Galactic Center region - down to ^1 mCrab 
level IjRevnivtsev et al. 2004|l . the Sagittarius Arm re- 
gion - down to ~1.5 mCrab level (Molkov et al. 2004) 
and the Galactic plane - down to ~2-3 mCrab level 
(|Bird et al. 2004|) were published. Here 1 mCrab corre- 
sponds approximately to ~ 1.4 x 10~ n erg s _1 cm -2 for 
a source with a Crab-like spectrum. Much INTEGRAL 
data has became public and this opens the possibility to 
construct a highly sensitive flux limited sample of sources 
in the inner part of the Galaxy and try to make a system- 
atic study of their broadband X-ray properties and their 
distribution in the Galaxy. 

In this paper we focus on the analysis of the popula- 
tion of high mass X-ray binaries located in the Galactic 
plane from the Norma to Sagittarius spiral arms. This 
Galactic region was selected because it has the best statis- 
tics of the available INTEGRAL data. We plan to extend 
this study towards further Galactic spiral arms when more 
INTEGRAL data are available. We present the high qual- 
ity broadband (1-100 keV) spectra of a number of newly 
discovered absorbed objects and report the discovery of 
X-ray pulsations from IGR J16465-4507. 



2. Data reduction 

For the reconstruction of the source broadband spec- 
tra we used all available X-ray data, which include data 
of the INTEGRAL observatory (all public data for rev- 
olutions up to 118 and the General Program data of 
the Galactic Center observations, Obs.ID 0120213 and 
0220133) as well, data of RXTE, ASCA and XMM- 
Newton observatories. The data of the latter observatories 
were reduced using standard tools recommended by the 
Guest Observer Facilities (http://legacy.gsfc.nasa.gov). 
The data of all INTEGRAL/IBIS /ISGRI observations 
were processed with the method described by Revnivtsev 
et al. (2004). In order to reconstruct the source spectra 
from IBIS/ISGRI data we used a ratio of fluxes measured 
in different energy channels to the fluxes measured by the 
ISGRI detector from the Crab nebula in the same energy 
bands. Detailed analysis of Crab nebula observations sug- 
gests that with the approach and software employed, a 
conservative estimation of uncertainty in measurements 
of absolute fluxes from the sources is about 10% and the 
shape of the spectrum is about 5%. The last value was 
added to the following spectral analysis as a systematic 
uncertainty in each energy channel. 

RXTE spectra of IGR J16318-4848 and IGR J16358- 
4726 were taken from the work of Revnivtsev (2003). In 
this work the special treatment was done in order to get 
rid of a strong contribution of the Galactic ridge emission. 

3. The sample and completeness 

The available INTEGRAL data now contain approxi- 
mately >20 Msec of observations, in particular ~>5 Msec 
of the inner parts of the Galaxy < 90°). In this paper 
we focus on the Galactic regions with Galactic latitudes 
|6| < 5° and Galactic longitudes from I > 325° to I < 50° 
(from the Norma spiral arm tangent to the Sagittarius spi- 
ral arm tangent). These data allowed us to obtain the time 
averaged map of the region and to construct a flux lim- 
ited sample of sources with a limiting sensitivity of ~1.5 
mCrab in the 20-60 keV energy band, that corresponds to 
an energy flux ~ 1.8 x 10 -11 ergs s" 1 cm~ 2 for a Crab-like 
spectrum. This criterion was chosen to provide a maximal 
survey area with best possible sensitivity. 

In the selected region of the Galactic plane (325° < 

1 < 50°) INTEGRAL detects 89 different sources with a 
time average flux of <^1.5 mCrab in the 20-60 keV en- 
ergy band. Among these sources there are 49 identified 
low mass X-ray binaries, 22 high mass X-ray binaries and 

2 active galactic nuclei. Also there are supernova remnants 
and single pulsars. Ten sources remain unidentified. A de- 
tailed analysis of the available INTEGRAL sky map will 
be presented in a separate paper. 

The list of high mass X-ray binaries presented in our 
sample is shown in Tabled In this Table we also present 
the value of the observed X-ray photoabsorption (column 
3), the value of the interstellar Galactic A% value inferred 
from radio maps QDickey fc Lockman 19 90), the average 
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Table 1. List of HMXB located in the Galactic plane from the Norma to Sagittarius spiral arms which were detected 
with INTEGRAL 



Source 


a (2000) 


5(2000) 


N H ,10 22 cm" 2 


N H ,10 22 cm" 2 


Flux, mCrab a 


u 

Comments 








Observed 


Galactic 






4U 1538-522 


235.584 


-52.376 


1.6[1] 


0.96 


16.5+0.2 


P 


AXJ161929-4945 


244.871 


-49.758 


H 


2.19 


2.0+0.2 




IGR J16318-4848 


247.953 


-48.801 


310, 


2.07 


21.8+0.2 




AX J1631.9-4752 d 


248.009 


-47.859 


18, 


2.13 


13.0+0.2 


P 


IGR J16358-4726 


248.990 


-47.407 


40, 


2.20 


3.11+0.21 


P 


AX J 163904-4642 


249.757 


-46.676 


58, 


2.18 


4.63+0.21 


P 


IGR J16465-4507 


251.648 


-45.118 


72, 


2.12 


8.8+0.9 


P 


IGR J16479-4514 e 


252.032 


-45.206 


12, 


2.14 


3.22+0.20 




OAO 1657-415 


255.199 


-41.653 


40 [2] 


1.76 


68.2+0.2 


P 


4U 1700-377 


255.982 


-37.841 


2-100(3] 


0.74 


185.0+0.2 




EXO 1722-363 


261.286 


-36.280 


50 [4] 


1.50 


7.18+0.13 


P 


IGR/XTE J17391-3021 


264.802 


-30.329 


5-6 [5] 


1.37 


1.50+0.11 




AX J1749.2-2725 


267.335 


-27.511 


10[6] 


1.62 


1.51+0.11 


P 


IGR/SAX J18027-2016 


270.677 


-20.278 


1-1.5[7] 


1.04 


4.06+0.14 


P 


AX J1820.5-1434 


275.131 


-14.553 


13[8] 


1.65 


2.94+0.20 


P 


AX J1838.0-0655 


279.523 


-6.921 


9, 


1.86 


2.36+0.24 




GS 1843+00 


281.412 


0.891 


2.3[9] 


1.01 


4.49+0.20 


P 


XTE J1855-026 


283.873 


-2.597 


15[10] 


0.73 


11.4+0.2 


P 


4U 1901+03 


285.914 


3.215 


7 f 


1.03 


75.5+0.2 


P 


4U 1907+097 


287.401 


9.843 


3-8 [11] 


1.75 


13.2+0.2 


P 


X1908+075 


287.699 


7.598 


10-50(12] 


1.48 


13.4+0.2 


P 


SS 433 


287.950 


4.990 


_g 


0.76 


14.2+0.2 


BH 


XTE J1858+034 


284.693 


3.429 


6 [13] 


1.89 


15+1 


P 



a - in the 20-60 keV energy band; 

6 - X-ray pulsar (P), black-hole candidate (BH); 

c - phot oabsorpt ion column obtained from our analysis of the broadband spectrum (see Table 2) 
d - the same as IGR J16320-4751; 

e - source localization was improved in our analysis of the archival data of the ASCA observatory: RA = 16 48 m 07 s , 

Dec = -45° 12' 21" (an accuracy of 0.9' ); 

f - determined from the analysis of the public RXTE data; 

3 - central X-ray source is completely absorbed; the observed X-ray emission emerges from a hot outflowing jet 

fe.g lWatson et al. 1986t; 

References: [1] -IRobba et al. 200l1 [21 -[Chakrabarty et al. 2002| [3 ] - |Haberl fc Day 1992| [4] - ILutovinov et al. 2004al [5] 
- ISmith et al. 19981 [6] - ITorii et al. 19981 [7| -|Augello et al. 2003| [8] -|Kinugasa et al. 19981 [9] - IPiraino et al. 20001 [10] 
- ICorbet et al. 19991 [11] - IRoberts et al. 20011 1 121 - ILevine et al. 20041 1 131 - IPaul fc Rao 19981 



source flux in mCrab values and comments about the 
source nature. 

Most of the presented HMXBs (16 of 23) are accretion- 
powered X-ray pulsars, one is a black hole candidate, 
4U1700-377 is a well-known high mass X-ray binary with- 
out detectable pulsations. Other sources were included in 
the list because of the strong intrinsic absorption in their 
X-ray spectra that is not typical for low mass X-ray bi- 
naries (IGR/XTE J17391-3021 was included based on the 
identification of the compaion star as a late O-type star, 
ISmith 2004|> . One of the sources, IGR J16465-4507, was 
not detected in the INTEGRAL time averaged map, as 
it was discovered only in September 2004 (revolution 232, 
ILutovinov et al. 2004bfl . Follow-up observations with the 
XMM-Newton observatory revealed a strong absorption 
in the source spectrum (see Table [TJ. Moreover, the im- 



provement of the source position ( |Zurita fc Walter7"2 004 ) 
allowed observers to propose a blue supergiant star as 
its possible counterpart in the system ((Smith 2004|) . The 
timing analysis of XMM-Newton data allowed us to de- 
tect pulsations of X-ray flux from this source (see be- 
low). Therefore we included IGR J 16465-4507 in the list 
of HMXBs and have a sample of 23 sources in total. 

We can estimate the number of extragalactic (mostly 
AGN) sources we should detect in the studied region 
with our sensitivity limit. For this purpose we can use 
either the flux-number function of Krivonos et al. (2005) 
in the energy band 20-60 keV or use the more conven- 
tional flux-number function in the 2-10 keV energy band 
(e.g. lUeda et al.200! ISazonov fc Revnivtsev 2004|) with 
the rescaled sensitivity limit (assuming, for example, a 
photon index of AGN spectra of T ~ 1.7). In both cases 
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Table 2. Best-fit spectra parameters of 8 highly absorbed sources in the 1-100 keV energy band 



Source 


r 


Nh,10^ cm" 2 


E c , kcV 


xV(d.o.f) 


AX J161929-4945 


1.3±1.1 


14±4 


31±8 


0.64 


IGR J16318-4848 


1.0±0.5 


310±70 


22±1 


0.95 


AX J1631. 9-4752 


0.7±0.2 


18±2 


13±1 


1.17 


IGR J16358-4726 


0.7±0.5 


40±10 


16±5 


1.16 


AX J163904-4642 


1.3±1.0 


58±11 


11±1 


0.57 


IGR J16465-4507 


1.0±0.5 


72±6 


30 


0.68 


IGR J16479-4514 


1.4±0.8 


12±4 


32±2 


0.97 


AX J1838.0-0655 


1.5±0.5 


9±3 


> 50 


0.82 



the estimated number of AGN in our region is of the order 
of 3-4 sources. Two AGN have already been identified - 
GRS 1734-292 flPavlinsky et al. 19920 and IGR18027-1455 
IjMasetti et al. 2004|) . 



4. Absorbed sources 

The study of absorbed X-ray sources at- 
tracted much of attention with the launch of 
INTEGRAL <|Matt fc Guainazzi, 20031 IRevnivtsev 20031 
IRodriguez et al. 2003| IWalter et al. 20031 Foschini et al. 
2004. ICombi et al. 20041 etc.). INTEGRAL is able to 
detect and measure the spectra of sources which are 
practically invisible for low energy X-ray instruments. 
The population of sources that demonstrates a high value 
of photoabsorption is not so small as to be disregarded - 
about 19 sources out of 23, presented in Table 1 have a 
significant (> 5 x 10 22 cm -2 ) intrinsic photoabsorption. 
In Fig. n we present the distribution of the observed JVh 
values in spectra of sources in our sample. 

4.1. Spectra 

To reconstruct the broadband spectra of absorbed sources 
we combine all available data from ASCA, RXTE, XMM- 
Newton and INTEGRAL observatories. The energy spec- 
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Fig. 1. Distribution of the absorption column value ob- 
served in X-ray spectra of HMXBs of our sample. The 
dashed line denotes the area where systematic uncer- 
tainties caused by high interstellar absorption in some 
Galactic regions prevents accurate measurements the 
source intrinsic photoabsorption value. 



tra of 8 sources in the 1-100 keV energy band are pre- 
sented in Fig. [21 These sources are the least studied in the 
broad energy band binaries located in the inner part of the 
Galaxy; most of them were discovered or rediscovered with 
INTEGRAL in hard (> 20 keV) X-rays. Spectral points 
in a soft energy band (< 20 keV) obtained with non- 
simultaneous observations of ASCA, XMM-Newton and 
RXTE were renormalized to match the INTEGRAL/IBIS 
spectral normalization. Here we assumed that the shape of 
the source spectra is not variable and used time averaged 
data. 

The obtained spectra can be well described by a sim- 
ple powerlaw model modified by the cutoff at high energies 
and photoabsorption at soft X-rays that is typical for ac- 
creting neutron stars (see e.g. the spectrum of IGR/XTE 
J17391-3021 o bserved with INTEGRA L in the Galactic 
Center region. iLutovinov et al. 2005 b). Best fit parame- 
ters are summarized in Table Best fit models are shown 
in Fig. by solid lines. 

Note that for the source IGR J16465-4507 the spec- 
trum is based on the data of XMM-Newton observa- 
tions performed on Sept. 14, 2004 and the measurement of 
INTEGRAL (|Lutovinov et al. 2004b|l . The solid line rep- 
resents the spectral form typical of new absorbed sources 
- a power law with an exponential cutoff (this value was 
fixed at the energy E c — 30 keV to better approximate 
the INTEGRAL point) and modified with the photoab- 
sorption, determined with the help of XMM-Newton data. 

It is seen that most of the sources have a somewhat 
similar spectral shape and parameters, in particular the 
folding (cutoff) energy. Only one source (AX J18380-0655) 
demonstrates a high energy cutoff value in its spectrum 
above ~ 30 keV. Such high values of the cutoff are not 
typical for neutron star systems, and therefore this sources 
is likely a black hole candidate. All other systems demon- 
strate X-ray spectra typical of accreting neutron stars. 

4.2. Pulsations in IGR J16358-4726 and IGR 
J16465-4507 

For some of the newly discovered sources coherent pulsa- 
tions of the X-ray flux were detected - IGR J 16358-4726 
l|Patel et al. 2004|), AX J163904-4642 IjWalter 2004jl . AX 
J1631.9-4752 l|Lutovinov et al. 2005aj) . 
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Fig. 2. Broadband energy spectra of 8 highly absorbed sources. Data in the hard energy part of spectra (> 20 
keV) represent the INTEGRAL observatory measurements; data in standard X-rays were obtained by ASCA for AX 
J161929-4945, AX J1631. 9-4752, AX J163904-4642, IGR J16479-4514, AX J1838.0-0655, RXTE for IGR J16318-4848, 
IGR J16358-4726 and XMM-Newton for IGR J16465-4507. Best fit models are shown by solid lines. Residuals arc 
presented in the bottom panels. 
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Fig. 3. Folded light curve of IGR J16358-4726 obtained 
with INTEGRAL/IBIS/ISGRI in the 18-60 keV energy 
band 

The available data of INTEGRAL and XMM-Newton 
observatories allow us to confirm pulsations in the flux of 
IGR J16358-4726 and detect pulsations in the flux of IGR 
J16465-4507. 

INTEGRAL observed IGR J16358-4726 during its 
outburst in spring of 2003 HRevnivtsev et al. "20 03b 
IRevnivtsev et al. 2003cf) . The source lightcurve in the I8- 
60 keV energy band demonstrates clear pulsc-like varia- 
tions. The total exposure time of INTEGRAL observa- 
tions used to construct the light curve of IGR J16358- 
4726 at the peak of its outburst was about ~ 90 ksec. 
The best period, obtained with the epoch-folding tech- 
nique, is P — 5980 ± 22 sec (the error represents the 
90% confidence interval; it was determined from the boot- 
strap Monte-Carlo simulations). The pulse profile of IGR 
J16358-4726 is presented in Fig.|31 As in the softer ener- 
gies (see Pa tel et al. 2004|l the source pulse profile consists 
of one peak. The pulse fraction in the 18-60 keV energy 
band (which is defined as P = (I max —Imin)/ (Imax+Imin), 
where I max and I m in are intensities at the maximum and 
minimum of the pulse profile) is equal to P — 70±10%, the 
same as in the CHANDRA energy band. The detailed tim- 
ing analysis of the data of the XMM-Newton observations 
of IGR J 16465-4507 allowed us to increase the number of 
known accreting X-ray pulsars in the region of the Norma 
spiral arm. 

The light curve of IGR J16465-4507 obtained with the 
XMM-Newton observatory during TOO observations of 
the source on Sept. 14, 2004 is shown in Fig. Strong 
variations are presented in the light curve. The best fit 
period determined with the epoch-folding technique is 
P = 228 ± 6 sec; the pulse fraction was estimated as 
~ 50 — 60%. The folded light curve is shown in Fig. 01 
(lower panel). Note that zero flux bins in the original light 
curve are caused by the Poissonian statistics. The mean 
count rate during the pulsation minimum is ~ 0.03 cnts/s, 
or ~ 1 cnts per bin (40 sec). We could not search for pulsa- 
tions in the INTEGRAL data because they are currently 
not publically available. 
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Fig. 4. (upper panel) Part of the light curve of IGR 
J16465-4507 obtained by XMM-Newton (the time bin is 
40 sec), (bottom panel) Corresponding pulse profile of IGR 
J16465-4507 folded with the period of 228 sec. 

For other sources no pulsations have been detected yet. 
However, similarities of their spectral shape to that of con- 
firmed high mass X-ray binaries and the presence of the 
strong absorption in their spectra allow us to tentatively 
suggest that all they are HMXBs. 

5. Spatial distribution of the X-ray binaries 

The obtained sample of Galactic X-ray binaries allows us 
to study their spatial distribution in the comparison with 
the previous work of Grimm et al. (2002). In that work 
authors used the data of RXTE/ASM (energy band 1-12 
keV) and that less sensitive to sources with high intrinsic 
photoabsorption. In addition to that, our limiting sensi- 
tivity is slightly better than that of Grimm et al. (2002). 
Assuming a spectrum of a typical accreting X-ray pulsar 
with zero intrinsic photoabsorption (power law with pho- 
ton index r = 1.3 with a high energy cutoff at 30 keV) 
our sensitivity of ~ 1.8 x 10 -11 ergs s _1 cm -2 in the 20- 
60 keV energy band corresponds to ~ 1 — 2 x 10~ n ergs 
s" 1 cm" 2 in the RXTE/ASM energy band 1-12 keV (cf. 
~6x 10 -11 ergs s -1 cm~ 2 in Grimm et al. 2002). 

In the region of our study (325° < I < 50°) we have 
detected in total 23 high mass X-ray binaries and HMXB 
candidates. This is more than 2 times larger than the sam- 
ple used by Grimm et al. (2002) in the same region. In 
the region / ~ 340° our sample contains 6 sources, while 
Grimm et al. (2002) considered only one source. This is 
not a surprise because a significant part of our sample 
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Fig. 5. Number- flux function of the Galactic HMXBs de- 
tected by INTEGRAL /IBIS /ISGRI in the inner part of 
the Galaxy. The dashed line represents the shape of the 
HMXB number-flux function obtained by Grimm et al. 
(2002) 

consists of new sources discovered by INTEGRAL which 
have strong photoabsorption. Note however, that the to- 
tal sample of HMXBs used by Grimm et al.(2002) includes 
more than 50 HMXBs and is larger than ours. 

The number-flux distribution of HMXB detected with 
INTEGRAL is presented in Fig. |5] The dashed line repre- 
sents the number-flux relation obtained by Grimm et al. 
(2002). 

The angular distribution of high mass X-ray binaries 
along the Galactic plane is presented by a histogram in 
Fig. in the comparison with the distribution of LMXBs. 
To estimate the significance of the difference of the ob- 
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Fig. 6. Angular distribution of high mass X-ray binaries 
(HMXBs) and low mass X-ray binaries (LMXBs) in the 
inner Galaxy. The number of LMXB is devided by 2. 




Fig. 7. Face-on view of the Galaxy with overlayed den- 
sities of HMXBs shown in gray scale (from 1 to 6, see 
Fig. IHJ). The distance of 8.5 kpc of the Galactic Center 
from the Sun is assumed. The star denotes the position of 
the Galactic Center. The open circle with the dot denotes 
the position of the Sun. Solid lines show the position of 
the spiral arms by Taylor & Cordes (1993), the dashed 
line shows logarithmic spiral arms by Vallee (1995), see 
also Russeil (2003) 

served HMXB distribution from the LMXB one we used a 
Kolmogorov-Smirnov test in two ways: for the whole sam- 
ple (325° < I < 50°) and for the first (0° < I < 50°) and 
fourth (325° < I < 360°) Galactic quadrants separately. 
We made this division in order to put the regions of in- 
terest in the middle of trial interval, where the K-S test is 
more sensitive. We build cumulative distributions for each 
of these cases and found that the probabilities that the 
HMXB distribution differs from LMXB one are ~ 75% 
for the whole sample and ~ 96% and ~ 97% for first 
and fourth quadrants, respectively (the significance of the 
central peak for the LMXB distribution is much higher, 
> 99.9%). It is known (Grimm et al. 2002) that low mass 
X-ray binaries trace the stellar mass of the Galaxy and 
are concentrated in the Galactic center. From Fig. Eland 
the K-S test it follows that the high mass X-ray binary 
distribution is different to the LMXB one and not concen- 
trated at the Galactic center. HMXBs, being the young 
population of X-ray sources in the Galaxy should be re- 
lated to star formation regions (e.g. Galactic spiral arms) 
rather than to the region with high stellar mass density, 
and such a connection was indeed observed (e.g. Grimm 
et al. 2002, Grimm et al. 2003). 

In Fig. we visualize the distribution of HMXB in the 
Galaxy with respect to the Galactic spiral structure. The 
number density of the sources is shown in gray scale. The 
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spiral model of the Galaxy is based on optical and radio 
observations of HII regions ( |Taylor fc Cordes 1993| ). The 
dotted circles denote arcs of circles centered on the Sun 
and show maximal distances at which our sensitivity limit 
(1.5 mCrab) allows us to detect X-ray sources with the 
marked luminosities (L x = 10 35 erg/s and L x = 10 35 - 5 
erg/s). The dashed lines denotes four logarithm spirals 
with pitch angle 12° adopted by Vallee (1995). This model 
predicts the position of spiral arms in the opposite (with 
respect to the Galactic Center) part of the Galaxy where 
measurements in the optical and radio bands are scarce, 
but these parts can be visible in X-rays. 

The detailed comparison of the density of the high 
mass X-ray binary distribution with the spiral structure 
of the Galaxy (regions of the enhanced star formation) 
is not possible at the moment due to the limiting statis- 
tics of the sample used. This is the subject of a sepa- 
rate paper in which more INTEGRAL data for the whole 
Galaxy will be used. Here we only mention several possi- 
ble reasons that may affect the observed HMXB density. 
They are: unknown distances to the sources and, as a se- 
quence, their unknown exact positions; dynamics of the 
Galaxy and complicated motions of the Galactic spiral 
structure and stars from the moment of star formation to 
the X-ray source stage (e.g. a significant displacement of 
HMXBs with respect to the current position of the spiral 
arms was observed by CHANDRA in spiral galaxy M83 
ISoria fc Wu 2 003) ; the possible influence of previously un- 
seen parts of the spiral arms, etc. 

6. Unidentified sources 

In order to understand the possible influence of unidenti- 
fied sources in our sample on the obtained distribution of 
HMXB we construct their distribution along the Galactic 
plane with the same coordinate binning as in Fig. H3 

This distribution (Fig. |HJ) shows that the unidentified 
sources cannot affect our conclusion about the distribu- 
tion of high mass X-ray binaries in the inner part of the 
Galaxy. There is a small peak near the Galactic Center 
that may indicate that several of the unidentified sources 
are LMXBs. One or two of the unidentified sources may 
be AGN. 



7. Summary 

— We presented a sample of high mass X-ray binaries in 
the inner part of the Galaxy from the Norma spiral arm 
region to the Sagittarius spiral arm region (—35° < 
I < 50°). This sample is significantly larger than the 
HMXBs sample in the same region of the sky, based on 
RXTE/ASM data, used by Grimm et al. (2002). The 
main reason for this is the hard X-ray energy band 
(18-60 keV) used in the present observations which 
helped to reveal a considerable population of absorbed 
HMXBs. 

— The majority of HMXB sources in our sample are sig- 
nificantly photoabsorbed. 



40 



20 
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340 



Fig. 8. The angular distribution of unidentified sources of 
our source sample along the Galactic plane. 



— We performed for the first time a spectral analysis for 8 
strongly absorbed sources of our sample. Based on the 
spectral approximations we argue that the majority of 
these sources should be neutron star binaries with a 
high mass companion. The exception is AX J 183800- 
0655, which likely contains a black hole. 

— We detected pulsations in the hard X-ray flux (I8- 
60 keV) from ICR J16358-4726 with the period P = 
5980 ± 22 sec, thereby confirming the results of the 
CHANDRA observatory IjPatel et al. 2004jl . In addi- 
tion we found pulsations with the period P = 228 ± 6 
sec in the X-ray flux of ICR J16465-4507 using the 
data of XMM-Newton. 

— We constructed the angular distribution of high mass 
X-ray binaries in the inner part of the Galaxy and 
showed that this distribution differs from the distribu- 
tion of LMXBs, which are concentrated in the Galactic 
Center. 

— A large fraction of the accretion powered X-ray pulsars 
in the Be-systems are transients. Therefore the con- 
tinuation of INTEGRAL observations of the Galactic 
plane might lead to the discovery of many more similar 
systems and considerably increase the statistical signif- 
icance of the observed displacement of the HMXB pop- 
ulation with respect to the spiral arm tangents. Data of 
CHANDRA and XMM-Newton on other galaxies and 
detailed study of systems in our Galaxy might provide 
additional information about the motion of stars and 
spiral waves. 
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